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Abstract 

The pion masses and the pion decay constant are calculated to 1-loop order in SU(2) twisted mass 
Wilson chiral perturbation theory, assuming a large pion mass splitting and tuning to maximal 
twist. Taking the large mass splitting at leading order in the chiral expansion leads to significant 
modifications in the chiral logarithms. For example, the result for the charged pion mass contains 
a chiral logarithm that involves the neutral pion mass instead of the charged one. Similar modifica- 
tions appear in the results for the neutral pion mass and the decay constant. These new results are 
used in fits to lattice data obtained recently by the European twisted mass collaboration. The data 
can be fitted well, in general better than with the standard chiral perturbation theory expressions 
that ignore the mass splitting. The impact on the extraction of low-energy couplings is briefiy 
discussed. 

PACS numbers: ll.15.Ha, 12.39.Fe, 12.38. Gc 



1 



I. INTRODUCTION 



Lattice QCD with twisted mass Wilson fermions jl|, Q has some advantages compared to 
its counterpart with a standard mass term. The most prominent one certainly is automatic 
0(a) improvement at maximal twist jsl. An efficient algorithm 0] allows simulations with 
sufficiently small pion masses to be in the chiral regime of QCD. Many results have been 
obtained in the quenched approximation as well as for unquenched 2-flavor-QCD (for a 
review see ref. j^). Recently, ffist results with additional active strange and charm quarks 
(2+1+1 flavor simulations [gI) have been reported too 

A drawback of the twisted mass formulation is the breaking of isospin symmetry, most 
clearly seen in a mass splitting between the charged and neutral pions. This breaking is not 
a fundamental concern; it is a lattice artifact of O(a^) and isospin symmetry is automatically 
restored in the continuum limit. Nevertheless, at nonzero lattice spacing the mass splitting 
is rather large. For example, for the 2+1+1 flavor simulations at a ~ 0.078fm one finds 
M^o/M^± ^ 0.54 for M^± ^ 320MeV and M^o/M^± ^ 0.77 for M^± ^ 400MeV.^ The 
neutral pion masses have a ten percent error; still, the neutral pion mass is roughly half as 
heavy as the charged one for the smaller charged pion mass. This is not a small effect. 

A large splitting in the pion masses is worrisome because of the chiral extrapolation 
which needs to be performed in order to get results at the physical pion mass. The standard 
tool here is chiral perturbation theory (ChPT) j8-10|. This low-energy effective field theory 



of QCD provides the quark mass dependence of observables. In particular, it predicts a 
characteristic non-analytic quark mass dependence, so-called chiral logarithms (chiral logs 
for short). In continuum ChPT without isospin breaking the chiral logs involve degener- 
ate charged and neutral pion masses. The same is true for the Wilson ChPT (WChPT) 



expressions in ref. HI], which incorporate the lattice spacing corrections assuming to be 
in the GSM regime where /i ~ qAqqi). In fact, at maximal twist the NLO expressions in 
the GSM regime are identical to those in the continuum. However, if the mass splitting is 
large, one can expect results involving the logs M^±lnM^±/A^ and M^olnM^o/A^. This 
modification may have a non-negligible impact on the chiral extrapolation. The extraction 
of the Gasser-Leutwyler (GL) coefficients, which are associated with the chiral logs, may be 
strongly affected by a large mass splitting. It is even possible that the modifications are so 
severe that the lattice data are not described at all by the continuum ChPT results. 

All this is reminiscent of staggered fermions with taste symmetry breaking. The lattice 
simulations by the MILC collaboration with Asqtad staggered fermions show a sizeable mass 
splitting between the various taste partners.^ For example, the ratio of the heaviest (taste 
singlet) pion mass M^r^ and the lightest (Goldstone) pion mass M^^ is approximately 0.77 at 
a ~ 0.09fm for a Goldstone pion mass M^, ^ 320MeV jl3|.^ These large taste splittings have 
been included in the chiral logarithms |lj4l6| (so-called staggered ChPT), and the lattice 
data clearly favor these modified logs. In fact, the continuum ChPT expressions cannot be 



^ Numbers from tables 5 and 8 of ref. M. 

^ For a comprehensive review see ref. |12| 

Numbers from tables III and IV of ref. hj. 
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fitted at all to the lattice data. 

In this paper we compute the pion masses (both for the charged and the neutral pion) 
and the decay constant to 1-loop order in SU(2) WChPT with the mass splitting taken into 
account in the chiral logs. In the language of WChPT we work in the Aoki regime 17| (also 
called LCE regime). For simplicity we work at maximal twist only, which is the relevant 
case in practice. As expected, we do find deviations from the continuum ChPT results. 
For example, the 1-loop result for the charged pion mass has a chiral log proportional to 
M^o InM^o/A^, i.e. it involves the neutral pion mass. The continuum chiral log proportional 
to M^± lnM^±/A^, on the other hand, is not present at all.'^ Similar changes are found for 
the neutral pion mass and the pion decay constant. 

There is a second aspect to these modifications that is equally important. A smaller neu- 
tral pion mass implies larger finite volume (FV) corrections. Typically, the FV corrections 
are exponentially suppressed with M^-L, where L is the spatial extent of the finite volume 
18N20|. If the neutral pion is significantly lighter than the charged one, the associated FV 
corrections are significantly less suppressed. Formulated the other way around one can say 
that the FV corrections due to the neutral pion mass are exponentially enhanced. This 



has been pointed out recently in ref. 2l|, where the enhanced FV corrections have been 



computed using the resummed Liischer formula |22|-|24l| . This enhancement can be large. 
For example, consider again the aforementioned lattice data with M^o/M^± ^ 0.54. The 
volume is such that M^±L ^ 4.0, which implies M^oL 2.2 and significantly larger FV 
corrections. 

In the last section of this paper we show that the most recent lattice data of the ETM 
collaboration are indeed sensitive to the modifications in the chiral logs. We reanalyze the 
data in ref. [3] using our new results. Indeed, the data prefer the results with a neutral pion 
mass smaller than the charged one, and the pion mass splitting is compatible with the one 
directly measured. More importantly, the central values for some of the extracted LECs 
change sizably with our new fit formulae. 

The rest of this paper is organized as follows. In section II we briefly summarize various 
results of twisted mass WChPT that we need later on. The primary purpose is to settle 
our notation and to prepare the 1-loop calculation presented in section III. The following 
section IV contains the analysis of the 2-1-1+1 flavor data in ref. 0] using our newly derived 
results. Final conclusions are drawn in section V. 



II. TWISTED MASS WILSON CHPT 



Wilson ChPT (WChPT) |25|, l26| is the low-energy effective theory for lattice QCD with 
Wilson quarks. It is based on a joint expansion in small pion momenta, quark masses and 



small lattice spacings a. Pedagogical introductions to WChPT are given in refs. 27|, l28 



Also, this finding has its analogue in staggered ChPT. The 1-loop result for the Goldstone pion mass has 



chiral l ogs involving the taste partners M- 
absent 151. 



The naively expected chiral log involving M^^ is 
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where many references to the original hterature can be found as well. 

WChPT has two sources of explicit chiral symmetry breaking, the quark mass and the 
lattice spacing, and the appropriate power counting depends on their relative size. The 
literature distinguishes two different regimes which seem to be the relevant ones for present- 



day lattice simulations. The GSM {generically small masses) regime UJ] assumes that the 



breaking of chiral symmetry due to the quark mass and the lattice spacing is of equal size. 



aAq^Q. The Aoki or LCE {large cut-off effects) regime 



29 



on the other hand, 
2 ^ 3 This is the 



assumes smaller quark masses or larger lattice spacings such that m ~ d^I^Q-^ 
regime we focus on in this paper, because it implies a pion mass splitting of the order of the 
charged pion mass itself. \a the following we collect a few results that have been published 
before in various papers 
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301 . 131j : the main purpose is to introduce our notation. 



The leading order (LO) chiral lagrangian in the LCE regime reads 



Here C2 denotes the standard LO lagrangian from continuum ChPT [ol. Eo|. 



£2 = ^(S.ss.st) 



(1) 



(2) 



(. . .) stands for the trace over the flavor indices. S denotes the usual SU(2)-valued Goldstone 
boson field which involves the pion fields in the standard way, specified explicitly below. / 
and B are the familiar LO low-energy coefficients (LECs).^ M denotes the quark mass 
matrix which in the presence of a twisted mass fi has the form 



M = m + ifias 



(3) 



where (T3 is the third Pauli matrix. Note that the standard (untwisted) mass m refers to the 
so-called shifted mass which includes the 0(a) shift to the additive mass renormalization 
25 1 . The lagrangian £^2 contains the leading O(a^) correction jssj and reads (we follow the 



notation of ref. ll| ) 



(4) 



Wqq is the LEG associated with the O(a^) 



a LEG of mass dimension three 26 



correction and a is defined by a = 2Woa. Wq is 
hence a has dimension two and W^g is dimensionless. 
The SU(2) fiavor symmetry is explicitly broken by the twisted mass term 'pointing' into 
the (T3 direction. Consequently, the ground state Evac, i-e. the minimum of the classical 



potential energy, is no longer equal to the identity but of the form 33M35 



exp^z^cTsJ 



(5) 



where is called the vacuum angle. It is determined by minimizing the potential energy 
which reads 



V 



lb 



ST) 



t\2 



(6) 



^ With our conventions the decay constant J^^KifK, 93MeV. 
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For notational convenience we introduced the short hand notation^ 



C2 



P 



(7) 



This combination of LECs is of mass dimension four. Usin g th e ansatz (jS]) in V we can 

(8) 



derive the gap equation (W/d(j) = 0, which can be written as 30 

2Bfi cos = sin (2Bm — 2c2a^ cos 



This equation determines the vacuum angle as a function of the variable parameters in the 
theory, the two masses m, fi and the lattice spacing a: (f> = (f>{m,^,a). Solutions of the 
gap equation and th e ^sociated phase diagram of the theory are discussed at length in 

It turns out that there are two qualitatively different scenarios 



the literature 25 



33|, |35| 



depending on the sign of C2. For C2 > there exists an Aoki phase 36| for fi = 0. Parity 
and flavor are spontaneously broken in this phase, and the charged pions are massless due 
to the spontaneous breaking of the flavor symmetry. Negative values of C2, on the other 
hand, imply a first order phase transition at m = 0. The pions are always massive except 
for the neutral pion which becomes massless at the end of the phase transition line, given 
by \2Bi^i\ = -2c2a^. 

As usual, the tree level pion masses are calculated by expanding the field S around the 
vacuum configuration. We parametrize S according to 



S(x) = SV2exp 5]z7r,(xK// 



(9) 



Using this form in ([6]) and making use of the gap equation the LO pion masses are found as 

iol 



Ml 



2Bfi 



ML 



ML 



AMI 



AMI 



2c2a^{l-t^ 



where we introduced 



t 



cos I 



(10) 
(11) 

(12) 



However, for some special values of t one has to be careful. Obviously, the charged pion 
mass seems ill-defined for t = ±1, but the gap equation immediately tells us that t = ±1 
is a solution only if /i = 0. This is the standard untwisted case and one finds (for positive 



M\ 



2Bm + 2c2a^ = ML 30|. Another case that 



C2 and outside the Aoki phase) 
requires care is t = and C2 < 0. It seems that the squared neutral pion mass can become 
negative for sufficiently small /i values. However, in this case the result (fTTl) is only valid for 



^ Note that our definition for C2 is not exactly the same as in [25|. It differs by a factor of pa^ and we 
liave dropped the terms proportional to the quark mass. 
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\2Bfi\ > — 2c2a^, keeping M^o always non-negative. As mentioned before, |2i?/i| = — 2c2Ci^ 
corresponds to the endpoints of the phase transition hne with vanishing neutral pion mass. 

So far the results in (fTOj) and (fTTj) are valid for arbitrary m, n and a (recall the implicit 
dependence on these parameters via t{m,fj,,a)). Usually one is interested in the special 
case of maximal twist only. Here maximal twist refers to the situation where the untwisted 
mass assumes a particular ('critical') value, m = rricr- The most popular and widely used 
definition specifies maximal twist as the point where the PCAC quark mass mpcAc vanishes 
(so-called PCAC mass definition). The PCAC quark mass is defined by 

{d,Al{x)P\y)) 
""^^^^ - 2(P^(x)P^(y)) ^^^^ 

with flavor index c = 1,2. The expressions for the axial vector current A^^ and the pseudo 



scalar density P'^ have been worked out in refs. [111 . l37l |. To LO these are the familiar 
expressions from continuum ChPT plus corrections proportional to powers of the lattice 
spacing.^ As any observable, the PCAC quark mass is a function of m, and a. Hence, cast 
into an equation the PCAC mass definition for maximal twist reads, reads 

mpcKcim = m^j., n, a) = 0. (14) 

Note that nicr does in general depend on the twisted mass and the lattice spacing: m^j.{^, a). 
The PCAC mass is easily computed at tree level [III, [sol, 31|- To LO in the LCE regime 



one finds 30 



r«pcAc = m- . (15) 

This vanishes for m = (which implies t = 0). For this particular value the result f lTOj) for 
the charged pion mass turns into the one from continuum ChPT, with m replaced by /i. 
The pion mass splitting between the neutral and charged pion mass is maximal in this case. 



Taking into account higher order terms in the calculation of the PCAC mass [11|, [SJJ, 138 
one finds that a vanishing PCAC mass implies t = 0(a). Dropping all terms higher than 
linear in a we can approximately write t = Xa, where X is some combination of LECs. 
Note that this result still guarantees automatic 0(a) improvement. Using t = Xa in (ITOl) . 
(ITT]) and expanding the square root we find small O(a^) corrections to the result for t = 0. 
Moreover, the same is true for the terms stemming from the 0(ap^, a/x) terms in the chiral 



lagrangian [ll|, l3l|. All these terms are associated with one power of t. Hence, at maximal 
twist with t = 0(a) these terms become scaling violations of O(a^), in agreement with 
automatic 0(a) improvement at maximal twist. 

In the following we always assume to be in the LCE regime. The precise meaning of 
this assumption is that the contributions and AM^ to the neutral pion mass in (fTTj) 
are of the same order, ~ AM^. At maximal twist this is equivalent to /i ~ C2a^/B. 
Assuming that all dimensionful constants are roughly equal to some power of Aqcd we 
recover the condition we have given before: /i ~ a^Aq^j-,. 



''' Note that the currents and densities can be given either in the twisted or the physical basis, which are 
related by a field transformation. Here we always work in the twisted basis. 
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III. PION MASSES AND DECAY CONSTANT TO ONE LOOP 



A. Propagators, vertices and counterterms 

Our goal is to compute the 1-loop corrections to the tree level pion masses given in (fTOj) 
and (fTTI) . In order to simplify the calculation we restrict ourselves to maximal twist, which is 
the relevant case in practice. We keep the L^? lagrangian at LO, hence the 1-loop calculation 
takes into account the following terms: 



LO : v^. a? 
NLO : v^M, M'\ fa\ Ma', a 



2^2 /i#^2 ^4 (16) 



Here the restriction to maximal twist implies already some simplification: In general the 
chiral Lagrangian contains also terms of 0(p^a, Ma) [ll| as well as 0(a'^). However, these 
are proportional to the solution t of the gap equation. Hence, at maximal twist with t = 0(a) 
these terms are promoted to 0(p^a^, Ma^, a^) terms, and the list given above is meant to 
include these NLO corrections. 

For the 1-loop calculation we need the propagators and interaction vertices stemming from 
the LO lagrangian. The propagators involve the tree level pion masses given in flTU]) and 
(ITT]) . In the 1-loop correction we can set t = 0, since t = 0(a) leads to O(a^) corrections (to 
the chiral logs) which are beyond NLO. Therefore, the (Euclidean space-time) propagators 
in momentum space read 

raft 

^"'(P') = 2 ^ > M| = 25/i, a, 6 =1,2, (17) 

G'\P') = -r^2 ' Ml = 2Bf, + 2c,a' . (18) 

The interaction vertices are the terms involving more than two pion fields. The four-pion ver- 
tices from the kinetic part in £2 are the same as in untwisted continuum ChPT. Expanding 
to quartic order in the pion fields we find 

^p^4. = ^ {{d.rr ■ Tif - {d.Trfix^) . (19) 

Here we use the short hand notation d^n ■ tt = ^^^^vTcTTc, (c^^vt)^ = J^c'^f^'^cd^Hc and 
analogously for tt^.^ Similarly, expanding the mass term and the lagrangian £^2 to quartic 
order we find (vr^ = (vr^)^) 

1 



^M,4^ = - 2^ ("^ COS yu sin 0) vr^ , (20) 
C-c2a^,iT, = — ^2c2a^ cos^ (pTi^ - sin^ vrVg j . (21) 



We still keep the summation over the Lorentz index implicit. 
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It is more convenient to express the quark masses in (120|) by the tree level charged pion mass 
and the LO pion mass splitting 2c2a^. With the help of the gap equation we can rewrite 
the sum of (EOl) and (EU) as 



UiAn + C,,a^An = -^^1 + ^2c2aV - ^2c2a\l - t^) n^nl . (22) 

The vertices proportional to o?t^ lead, after contracting two of the four pion fields, to 
corrections proportional to a^t^M^ In M^/A^. For maximal twist with t = 0(a) these are 
corrections higher than NLO, hence they are ignored in the following. We can also ignore 
the three-pion vertices stemming from the mass term, 

'0-M,3n = — (/icos0 - msin0)7rV3 . (23) 

Making again use of the gap equation this can be brought into 

CM,3. = -^tVr~^7r'n,. (24) 

In order to form a loop diagram contributing to the self energy of the pions one needs 
two of these vertices. Therefore, with t = 0(a) this correction is at least proportional 
to (c2a^)^a^ ~ a^, which is a NNLO correction in the LCE regime and beyond the order 
considered here. 

We conclude that for an NLO calculation at maximal twist we can set t = in (!22|) . and 
there is only one extra vertex proportional to 2c2a^ besides the familiar ones from continuum 
ChPT. 

Finally, we need the analytic NLO corrections which also provide the necessary counter- 
terms for the divergent contributions of the loop corrections. The full NLO lagrangian in 
the LCE regime at maximal twist consists of the following parts: 

-^NLO = C4 + Cp2a'2 + CMa^ + + ^p'^a + ^Ma + ^a^^ ■ (25) 



The first part £4 denotes the standard NLO lagrangian of continuum ChPT [9[ (we again 



follow the notation of 11|). 

A = L45{df,Edf,E^){ME^ + SM^) - Lq^{MT) + SM^)^ , (26) 

with 

M = 2BM (27) 

and the mass matrix M defined in ([3]). Note that we dropped all terms in £4 that do not 
contribute to the pion masses (and decay constant), in particular the 0{p^) terms involving 
the GL coefficients Li,L2. 

The next three parts in £nlo are constructed in appendix A of j29j. Although a degen- 
erate untwisted quark mass matrix was assumed in this reference the generalization to the 
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case with a twisted mass is straightforward. We find 



+ a3a2(9^(S + St))(a^(S + St)), 

The coefficients undetermined LECs. 

The remaining corrections stem from the lagrangians 



(28) 
(29) 
(30) 



£ 2 



nven in ref. Ill], and the O(a^) correction 



W^45a(S + St)(a^Ea^Et), 



(31) 
(32) 

(33) 



given in ref. 29 



Ahhough in general of lower order in the chiral expansion than the ones 
discussed so far, these corrections are promoted to NLO terms at maximal twist. Let us 
demonstrate this for the contribution to the pion masses. Expanding (E + E^) into pion 
fields we obtain 



(34) 



where the ellipses stand for the terms with three or more pion fields. Important is the 
factor t, which is of 0(a) at maximal twist. Setting t = Xa with some constant X we find 
a tree level correction of C^a to the pion masses proportional to WqsXo^M^. This is a 
NLO correction in our counting and should be taken into account. Similar arguments can 
be made for the lagrangians £^2^ and £^3. 

It is straightforward to expand all terms in fl26l) - fl33|) to quadratic order in the pion 
fields. This leads to the analytic NLO corrections to the pion masses and wave function 
renormalization, which are of the generic form y47r^/2, i?7r|/2 and C(9^7r)^/2, with A,B,C 
being some combinations of the LECs appearing in (l26l) - (l33l) . For practical applications 
(fits to lattice data) there is no need to keep track of the individual LECs. However, one 
should make sure that the LEC combinations in the final results for observables are linearly 
independent, and this is the reason why we quoted all terms contributing to £nlo- 

Finally, the last NLO correction to the pion masses stems from the LO results. Setting 
t = Xa in ( fTOl) and (fTTl) we find a 0(M|.a^) correction to the charged pion mass and 
0(MQa^,a^) corrections to the neutral pion mass. These are NLO corrections and need to 
be taken into account as well. 



It is a matter of taste whether one uses a or a in these expressions, the difference is just a multipHcative 
constant 2Wo. However, the mass dimension of the LECs depends on this choice, because a and a have 
mass dimension —1 and 2, respectively. 

We dropped a term hnear in 713 as well, which does not play a role here. 
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B. Pion masses to NLO 



With the resuhs given in the previous section the NLO calculation of the pion masses 
is straightforward. The 1-loop diagrams contributing to the self energy of the pions are all 
tadpole diagrams and stem from f|T9l) and fl22l) . Except for the vertex from the last term in 
fl22|) these vertices are just the ones that also contribute in a continuum ChPT calculation. 
The only (but crucial) difference is that one has to keep track of the flavor indices for 
the pion in the loop, because the charged and the neutral pion have different masses. All 
loop diagrams lead to the familiar scalar integrals, which are conveniently regularized by 
dimensional regularization. The divergencies are removed by introducing renormalized LECs 
at NLO.^^ For the charged pion mass we find the result 

M.^.,NLO = Ml (^1 + ^ In ^ - (l45 - 2L,s^ + C,a^^ . (35) 

M± and Mq are the LO masses in f|T7|) . f|T8l) and fi denotes the renormalization scale. Here the 
NLO LECs are renormahzed couplings and depend on that scale, L1^{fi) , Llg{fi) and C[(/i), 
but for brevity we drop the superscript and do not make the dependence on fl explicit. 

As already announced in the introduction, the chiral log in ( l35l) contains the neutral 
pion mass and not the charged one. If the mass splitting is large this is a non-negligible 
modification. However, note that we recover the correct continuum result if a goes to zero, 
since Mq — )■ M± in this limit. 

The LEG Ci is a combination of LECs and incorporates the 0(M|a^) corrections to the 
charged pion mass from various sources, as discussed at the end of the last section. It is 
easily checked that £^4 in (!30|) does not contribute to the charged pion mass, so there is no 
O(a^) shift in (!35|) . This is also expected qualitatively. Suppose C2 > 0. In this case there 
exists an Aoki phase at /i = 0, and the charged pion mass must vanish for /i — )■ 0. This 
excludes an O(a^) correction in the result for the charged pion mass. In fact, for the same 
reason there is no O(a^) shift in the LO mass in (fTOl) . 

In practical applications it is convenient to introduce A3, defined by 

16(L45-2L68) = ^ln||. (36) 
In terms of A3 we can rewrite (135|) according to 

M.^. NLO = Ml (1 + ^ In ^ + CM^a^) . (37) 



Besides the divergence proportional to e ^ the finite part ln47r — 7 + 1 is also subtracted. 
That the charged pion contains a chiral log involving the neutral pion mass has already been noted in 
ref. [3^. However, there is a discrepancy in the prefactor of the chiral log. The source of this discrepancy 
is currently looked for [4^. Note, however, that the result in ref. [s^ does not reproduce the correct 
continuum chiral log in the continuum limit. 
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The new combination of LECs, Cm±, involves Ci and L45 — 2^68- The scale dependence 
drops out in this combination and Cm± is scale independent. 

The calculation of the neutral pion mass is completely analogous, and we find 

M^o^^^o = Ml(^l + (^2Ml In ^ - In - ^M|(L45 - 2^68) + C^a 

+ 2c,a^l-^\n^ + C,a']. (38) 



This result contains both types of chiral logs, involving the charged and the neutral pion 
mass. Still, in the continuum limit a — >■ we recover the continuum result, as expected. 

C2 and C3 are combinations of LECs associated with the 0(M|.a^) and O(a^) correc- 
tions. As before, these renormalized coefficients are scale dependent and the superscript 
'V is suppressed. As for the charged pion mass, these coefficients represent various NLO 
corrections stemming from (l28l) - (!33|) as well as from the corrections to the LO pion mass 
in (fTTj) with t = 0(a). In contrast to the charged pion mass there is a correction of 0{a'^). 

It is again convenient to replace L45 — 2^68 by A3. Similarly, C2 and C3 can be traded 
for a dimensionful scale S3 and a dimensionless and scale invariant coefficient Cmq'- 

Mi,™-Mi(l + ^(2M|l„^-Mj,„|' 



+ 2c2a' (^1 - ^ + J • (^^) 

For completeness we also quote the NLO result for the mass splitting. Although one could 
directly take the difference of ( 139|) and (137|) it seems beneficial to start with the differences of 
(!38|) and (!35|) . The contribution proportional to L45 — 2^68 drops out and the two O(M^a^) 
contributions can be combined at this stage. Then one proceeds as for the neutral pion mass 
and finds 

, ^ Ml ^ ^ 



AM^Vlo = 2c2a^ ( 1 - ^ In ^ + Cka 



S3 and Ca differ from the analogous coefficients in f l39p by terms proportional to Ci. 



C. Decay constant to NLO 

An equally important observable besides the pion mass is the pion decay constant /tt, 
defined by the one-pion matrix element of the axial vector current. In this section we 
compute the decay constant /^ri of the charged pions to NLO. For simplicity we write 
= since we never encounter the decay constant of the neutral pion. We will find that 
the chiral log in is also modified by a large pion mass splitting. 
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In twisted mass QCD one usually does not compute by the matrix element involving 
the axial vector current. Instead, one makes use of the so-called indirect method [2I, |4l| where 
the decay constant is given by 

U = , (41) 

G. = |(0|P'^(0)|7r'^(p-))|, a = 1,2, (42) 
where P°'{x) is the pseudo scalar density. Relation fj^Tl) is based on the exact Ward identity 

d^V;^ = -2^e^"''P^ (43) 

involving the vector current V"^. At maximal twist the vector current corresponds to the 
physical axial vector current. This implies (HTjl . and its benefit is that one does not need to 
compute any renormalization factors like Zp or Za- 

The right hand side of (HTjl is straightforwardly calculated in WChPT. The pion mass is 
given in the last section, and the missing piece is the matrix element G^. The expression 



for the pseudo scalar density can be found in ref. 1]|. To LO it is the familiar expression 
from continuum ChPT, 

^Lo = ^(^a(S-St)). (44) 

Expanding to linear order in the pion fields we find the tree level result G7r,LO = fB. Using 
the tree level result (fTOl) for the charged pion mass we obtain 

fn,LO = fVr~^. (45) 

This is the result for the right hand side of (14T1) for an arbitrary twist angle. At maximal 
twist (t = 0) it turns into the well-known LO result for the decay constant. If t = 0(a) it 
receives, after expanding the square root, an O(a^) correction. Note that the origin of the 
factor a/1 — is the result for the pion mass. Gn^LO is exactly as in continuum ChPT and 
does not add any modifications due to the non-trivial ground state. 

In order to compute the decay constant to NLO we need the NLO expression for the 



pseudo scalar density in the LCE regime. Most terms can be taken from ref. where 
P^i^Q has been derived for the GSM regime. Missing are the O(a^) corrections, but these 
are easily constructed (see below). 

The general structure of at NLO is 

^NLO = ^Lo(l + A^NLO.GSM + AP,^). (46) 



The result correctly vanishes for t = 0. This corresponds to the untwisted case. The Ward identity (|43|) 
still holds (for /i = 0), however, is the physical vector current and its one-pion matrix element vanishes. 
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The correction APnlo.gsm can be found in eq. (47) of (where it is denoted by V) 



a4 



AP^,o = -^(5,S5,St) + ^(Mts + stM) + + St) . (47) 

These terms are the corrections of 0(p^,M, a) to the leading 1 in f H6|l . The correction 
proportional to L45 contributes at NNLO only and can be ignored in the following. In the 
mass term we can set the untwisted mass m to zero since we are interested in the case of 
maximal twist only. The correction proportional to aWes is effectively a O(a^) contribution 
at maximal twist because of the contribution (S + E^) (recall the discussion after eq. flM|) ). 



The method used in ref. ll| for the construction of the pseudo scalar density is easily 
extended to the O(a^) terms in APa-2. The pseudo scalar density is obtained by a functional 
derivative of the effective action with respect to the mass, which is promoted to a space-time 
dependent spurion field in an intermediate step. For example, the 0(a) correction in M7\\ 
is obtained from the Cmu lagrangian in (l32l) . In complete analogy the lagrangian Cmu'^ in 
( l29l) yields the O(a^) corrections 

AP,. = -— (6i + 62(S + St)2j. (48) 

The term proportional to 62 is effectively an O(a^) term because of the factor (S + S^)^. 

Dropping all terms that contribute beyond NLO only we can use the following (incom- 
plete) NLO expression for the pseudo scalar density: 

^NLO = Plo (^1 + ^M|(za3(St - S)) - j,h,osa'^ . (49) 

where we introduced the effective LEG bi^eS = bi — SWqsWoX which includes the remnant 
O(a^) correction from the 0(a) term in ( H71) (as before we have set t = Xa)}^ 

The 1-loop calculation of the matrix element is now standard. Expanding Plq in ( H9ll 
up to cubic order in pion fields one obtains the terms that lead to the 1-loop corrections. 
Since Plq is the familiar continuum expression, the calculation is as in continuum ChPT, 
except for the fact that one has to keep track of the flavor index of the pion in the loop since 
the masses are different. The result of the calculation is 

1 - + -^(4^68 - L,5) + C,a^ j . (50) 

We introduced C4 as a short hand notation for the contributing combination of LECs. 
It contains bi^ca and also ai of the lagrangian £a2p2 (it contributes to the wave function 



14 The LEG W in ref. 
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corresponds to Wes in our case since we directly started with the chiral lagrangian 



parameterized in terms of the shifted mass. 
15 Wc remark that the expression of the pseudo scalar density is determined by the Ward identity (j43|) . 
In lattice QCD many pseudo scalar densities can be defined, all differing by 0(a). The corresponding 
expressions in WChPT differ too, since they have to accommodate these differences. This is analogous to 



the vector and axial vector currents, which have been discussed in ref. 37 j 
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renormalization). Forming the ratio in (HT]) with the result (135|) for the charged pion mass 
we finally obtain the NLO result for the decay constant: 



( I f . Ml „ Ml\ ml , , ' 



(51) 



Both (130]) and (I3T]) turn into the known continuum ChPT result for a — )■ 0. Also, the 
GSM result at maximal twist [ll| is correctly reproduced. The analytic lattice spacing 
dependence is O(a^), as expected from automatic 0(a) improvement. However, the chiral 
log is not the naively expected one: the averaged chiral log \M\ In Ml + Mq In Mq ] / 2 appears 
in the result. 

As for the pion masses it is convenient to slightly rewrite the result ( 15T|) . In terms of the 
scale A4, defined by 

8^45 = ^ In ^, (52) 
the result for the decay constant turns into 

/.,NLO = / (^1 - (m| In ^ + Ml In ^ + C,a^ . (53) 

In analogy to Cm± we introduced the notation C/ for the final combination of O(a^) LECs. 

We followed the actual numerical computation of the decay constant and calculated /tt 
with the indirect method. Alternatively, one can also compute the standard matrix element 
involving the (physical) axial vector current. The result must be the same, of course. We 
partially confirmed this by calculating the 1-loop correction for the axial vector current 
matrix element. We indeed found the chiral logs as in ( ISTl) . This also provided a non-trivial 
check for the pion mass computation in section IIII B[ 



D. Finite volume corrections 

So far all the calculations were done in infinite volume, but the finite volume (FV) correc- 
tions are known 18N20| and easily included. Most relevant in practice is a finite space-time 
volume with geometry x T and periodic boundary conditions in each direction. As has 



been shown in ref. 20[, the chiral lagrangian is as in infinite volume and the finite volume 
dependence enters through the pion propagators only. The position space propagator, ob- 
tained by Fourier transforming ( IT71) and (ITSl) . involves a sum over the discrete pion momenta 
instead of an integration. 

As a concrete example we quote here the modifications of the previously derived results 
for a finite spatial volume L^, assuming the temporal extend T to be much larger so that it 



This is also reminiscent of the 1-loop result in staggered ChPT: The pion decay constant involves the 
average of all sixteen chiral logs that one can form with the available taste partners 16 1. 
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can be taken infinite. In this case the FV corrections are included by the simple replacement 

3^' 



, , M2 „ ^^^^^ 

^ \n — + g,{ML) 



(54) 



UML) 



4 
ML 



/Ci(|n|ML) 



\n\ 



+ 77-2 + n| , 



(55) 



in every chiral log. Ki is the modified Bessel function of the second kind. The sum runs over 
all triples n = (ni, 'n.2, ris) where the are integers. Taking into account the degeneracies in 
this three-dimensional sum it collapses to a simple one-dimensional sum j23, 42|. The Bessel 
function Ki is exponentially small for large arguments, Ki{ML) ^ exp(— ML) / y/ ML, hence 
the sum in (155!) converges rapidly and the FV correction ^i(ML) is small. 

The results of the previous section contain chiral logs with the charged and the neutral 
pion mass. We therefore obtain FV corrections governed by the arguments M±L and MqL. 
The latter are significantly (exponentially) larger if the neutral pion mass is much smaller 
than the charged one, as has been pointed out in ref. 
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IV. NUMERICAL ANALYSIS 

The calculations of the previous section were triggered by the numerical results of the 
ETM collaboration in ref. As already mentioned in the introduction, measurements 
of the charged and neutral pion masses showed that the latter is significantly lighter, with 
M^o/Mn-i approximately 0.54 and 0.77 for two of the simulated quark masses at a ~ 0.078fm. 
These values correspond to roughly 0.71 and 0.41 for the ratio 

^ ^ — ^ f« ' ^ ' , (56) 

and this tells us that at least these data points are in the LCE regime. It is therefore 
interesting to check whether the data show evidence for the modified chiral logs and, provided 
this is the case, whether the extraction of Gasser-Leutwyler coefficients is affected by using 
the results of the previous section instead of the continuum ChPT results. 

Data for the charged pion mass and the pion decay constant is given in table 8 of ref. 
[3]. Measurements have been made for two different lattice spacings (a ~ 0.08585fm and 
0.0782fm) and for various (6 and 5) different light quark masses. The charged pion mass 
ranges roughly between 270 and 510MeV. In total there exist 22 data points for a combined 
fit to the data. This is rather large compared to the number of fit parameters in the 
expressions fl37|) and fl53|) : Four continuum ChPT parameters /, i?,A3,A4, and three more 



14| denotes the FV correction by 6i. The notation gi is useiLin refs. 23l l42l|. 



1^ Ref. 

A similar analysis of quenched lattice data can be found in ref. [43 
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associated with the nonzero lattice spacing, C2,Cm± and Cf. Even if one drops a few data 
points for the heavier pion masses there are still enough data points to perform a fit. 

All fit results presented in this section were obtained by fits to the dimensionless lattice 
data for a/^r and the ratio 

fi=i^^. (57) 

where /io denotes the bare twisted quark mass.^^ The renormalized quark mass fi used in 
the previous section renormalizes multiplicatively, thus fi is proportional to /iq and we have 
Bfi = BofiQ. Instead of quoting the fit results for A3, A4 we give the results for 

^3,4 = ln(-^), (58) 

where phys = 139.6MeV denotes the physical pion mass. 

Results of fits to the data with the pion mass splitting ignored in the chiral logs have 
been reported in ref. 0]. We also give results for such fits, but the results are expected to 
differ slightly for two reasons. In ref. 0] the resummed FV corrections of 23 1 were used, 



which differ from the FV formulae of section UlI Dl even if we set C2 equal to zero.^*^ Another 
difference concerns the error analysis. The statistical errors for the fit parameters in ref. 0] 
are estimated by generating bootstrap samples from the bare correlation functions for the 
pion mass and decay constant. This cannot be done having only the data in table 8 of [3] 
available. Instead, the fit results and error estimates given here were obtained by fitting the 
data with the Levenberg-Marquardt algorithm 4^ . 

In order to check for these potential differences we performed separate fits to the data 
at fixed lattice spacing with Cm± and Cf set to zero. The results are given in table [D 
We performed three fits that differ in the data points included in the fit. Fit I included all 
data points while for fit II (III) the data points at the smallest (heaviest) quark mass were 
excluded. For all fit ranges two fits were done, one that includes C2 as a fit parameter (right 
subcolumn) and one without, setting C2 equal to zero (left subcolumn, indicated by a dash). 

Fit II with C2 = has been done by the ETM collaboration and their results in table 
9 of j3] should be compared with ours (for the readers convenience we have collected the 
relevant fit results in table [TTI- The data points for the smallest quark mass were dropped 
by the ETM collaboration because the data did not fully comply with the tuning condition 
mpcAc/yLto < 0.1 for maximal twist (see section 3.2 in ref. 01). 

The mean values for / and the LECs I3 4 agree to a very good degree. Our errors for 
I34 are somewhat larger, which is not unexpected. The error for the decay constant, on the 



The conventions in ref. [7| are such that / « 130MeV. Therefore, in the results of the previous section 
the replacement — >■ /^/2 in the prefactor of the chiral logs has to be made. 



23| contain two more NLO LECs, and 12- 



Note that the formulae of 

For fits at one lattice spacing the constants Bq (/) and Cmj- (Cf) would not be independent if the chiral 
log was absent. Even in the presence of the (small) chiral log these fit parameters arc not very well 
determined individually. 
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3 59l'1 1 3 11 ^94^1 




4.74(2) 4.71(4) 


4.74(3) 4.71(7) 


4.74(4) 4.70(6) 


-2c2a2 (MeV^) 


[214(27)]2 


[213(90)]2 


[215(29)]2 




12 12 


10 10 


10 10 


X^/^-dof 


8.81/8 5.48/7 


6.50/6 5.46/5 


8.36/6 5.34/5 




0.55 0.86 


0.59 0.71 


0.40 0.72 



Fit {(3 = 1.95) 


I 


II 


III 


Fit range: a/io,min 


0.0025 


0.0035 


0.0025 


'^/^0,max 


0.0085 


0.0085 


0.0075 


maximal M^± (MeV) 


496 


496 


461 


2Boa 


4.94(4) 4.72(8) 


4.96(6) 4.61(8) 


4.98(5) 4.83(18) 


f (MeV) 


119.9(1.3) 128.0(3.3) 


120.5(1.4) 134.3(3.5) 


119.3(1.3) 124.3(6.4) 




3.66(8) 3.36(18) 


3.70(10) 3.16(23) 


3.74(9) 3.60(24) 


k 


4.67(3) 4.62(5) 


4.65(3) 4.52(7) 


4.68(3) 4.66(6) 


-2c20? (MeV^) 


[190(31)]2 


[249(22)]2 


[150(89)]2 


'^data 


10 10 


8 8 


8 8 


X^/f^dof 


13.3/6 10.6/5 


11.2/4 6.3/3 


4.1/4 3.80/3 


Q 


0.10 0.23 


0.08 0.39 


0.67 0.70 



TABLE I: Results for the fits to the data at fixed lattice spacings, a ~ 0.086fm (/? = 1.9, top) 
and a ~ 0.078fm (/3 = 1.95, bottom). Right subcolmnns correspond to fits with ci as a free fit 
parameter, left subcolumns to fits with C2 set to zero. The parameters Cm±, Cf are always set to 
zero in these fits (see text). 

other hand, is slightly puzzling: Our errors for / are roughly at the one percent level, while 
the error in ref. |7| is more than an order of magnitude smaller. The reason for this strong 
discrepancy is not clear to us. 

Let us turn to the fits in table [T] that include C2 as a fit parameter. Qualitatively we can 
say that the data sets prefer a negative value for C2 with — 2c2a^ ~ (200MeV)^. Although 
the error is quite large the sign is in agreement with a neutral pion lighter than the charged 
ones. All fits with C2 as a fit parameter have slightly better values for the x^/''^dof and 
the goodness of the fit Q, but the improvement is not dramatic. Noteworthy is that the 
values for / are systematically larger, the ones for I3 systematically smaller compared to the 
fits with C2 set to zero. The worst discrepancy with a 3.9a difference appears in Fit II at 
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/3 = 1.90 


/3 = 1.95 


combined 


/ (MeV) 


120.956(70) 


121.144(83) 


121.031(54) 


h 


3.435(61) 


3.698(73) 


3.537(47) 


h 


4.773(21) 


4.673(25) 


4.735(17) 



TABLE II: Selected results of the fits performed by the ETM collaboration, taken from table 9 of 
Ref. The first two columns refer to the separate fits, the last one to the combined fit. 

P = 1.95 for /, but most differences are (roughly) between la and 2.5a. These differences 
are sizable and not negligible. However, before one can draw firm conclusions these results 
need to be corroborated by fits that properly take into account any correlations in the data. 

The main motivation for the WChPT calculations in the previous section is a combined fit 
to the data at both lattice spacings, since this amounts in a combined chiral and continuum 
extrapolation. The results of such fits are given in table lllli As before, various fits were 
done with different ranges for the pion masses. Fits I to III include the same data points 
as in the separate fits. Fits with even more data points excluded were done as well: The 
smallest and the largest pion mass (at each lattice spacing) were dropped in fit IV, while the 
largest two (three) pion mass data points were excluded in fit V (VI). All fits include Cm± 
and Cf as free fit parameters. The ratio of the two lattice spacings, on the other hand, is not 
a fit parameter but included as the fixed ratio = 0.0782/0.08585. Results in table UTTl that 
include the lattice spacing refer to the values at the larger lattice spacing a ~ 0.08585fm 
{(3 = 1.9). 

A combined fit requires the ratio of quark mass renormalization factors at the two 
lattice spacings. This ratio is not available to us so we set it to 1. We expect this to be a 
good approximation because the two lattice spacings are very close: The finer lattice spacing 
is less than 10 per cent smaller than the coarse one.^^ 

None of the results in table IIIII can be compared with results in ref. [3] • Although the 
ETM collaboration has done a combined fit, the coefficients Cm± and Cf were set to zero in 
this fit. 

From table IIIII we draw the following conclusions. 

1. Quite generally, all fits are satisfactory with respect to their and Q values, even 
the ones with all data points included. The largest values for x^/'"'dof is about 2 in 
fit II. The quality of the fits improve if data points are dropped. Comparisons of the 
fits I with II and III as well as II with IV shows that decreases substantially if the 
data for the heaviest pion mass are excluded. The improvement is less significant for 
dropping the smallest pion mass data. This is in agreement with ChPT as a low-energy 
effective theory. 



We also did fits with the renormalization faetor ratio ineluded as a free fit parameter. We obtained ratios 
between 0.98 and 1.02 with an error of about 0.02. The results for the other fit parameter agree with the 
ones given in table IIIII within the errors. However, the error estimates for the fit parameters 2_Boa and 
CM±a^ are larger by a factor 4 to 5. 
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11.1/8 9.70/7 
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0.52 0.64 


0.73 0.92 


0.93 0.95 



TABLE III: Results for the combined fits to the data for both lattice spacings. As in table HI right 
subcolumns correspond to fits with C2 as a free fit parameter, left subcolumns to fits with C2 set to 
zero. The values for the fit parameters involving the lattice spacing refer to a ~ 0.086 fm (/3 = 1.9). 



As already mentioned, the reason for dropping the data at the smallest pion mass by 
the ETM collaboration was the potential violation of tuning to maximal twist. It is 
not easy to exactly quantify a small mistuning, but the results in table IIIII show that 
including the data at the lightest pion mass is less influential than including the data 
at the heaviest pion mass. The central values for flts IV and V are in good agreement 
and either of it seems (at least to us) to be a good candidate for obtaining reliable flt 
results. 

2. Fits with flt parameter are always better than without. However, the im- 
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provement is less significant for the fits with the heavier mass data excluded. Never- 
theless, the data prefer a negative C2 in agreement with a neutral pion mass smaller 
than the charged one. The LO pion mass splitting — 2c2a^ is roughly (200MeV)^ (at 
a ~ 0.086fm) with a large error. Note that the errors given in table IIIII are slightly 
misleading: The square root of the fit parameter — 2c2a^ is quoted and its relative 
error is half as big. 

Note that the neutral pion mass M^o nlq is not predicted by the values in table IIIII 
because it depends on two extra parameters, H3 and Cmq- Hence, we cannot check 
whether our fit results are in agreement with the direct measurements of the neutral 
pion mass. Moreover, it is not possible to include the data for the neutral pion mass 
in the fit, since only two data points are given in ref. 0]. Once more data become 
available it will be very interesting to attempt combined fits with the neutral pion 
mass data included. 

3. Comparing the fit results with and without C2 as a fit parameter one can observe: 
(i) the central values agree within errors, (ii) the central values for / are somewhat 
larger, and smaller for the Gasser-Leutwyler coefficients 13^4, (iii) the errors are in 
general larger with C2 included. Whether these observations persist for smaller errors 
cannot be said here. An error analysis as the one in ref. [3] may lead to smaller errors 
and perhaps to a firmer conclusion. 

4. Our last observation concerns the analytic O(a^) corrections to the LO LECs, Bq and 
/. Consider, for example, fit II with C2 excluded from the fit. The results for Cm± 
and Cf in the combined fit mean that there is a 20% and 10% O(a^) error in Bo and 
/ at a f« 0.086fm. These numbers are consistent with the separate fit at this lattice 
spacing (fit II in top of table [T]). Note that we have set C^j^ and Cf equal to zero in 
the separate fit, so the correction associated with these LECs is effectively absorbed 
in Bq and /: / in the separate fits corresponds to /(I + a^C/) in the combined fit, and 
analogously for Bq. 

Our results here show that the separate fits can be quite misleading for the estimates 
in the continuum limit. For example, the central value for / is unchanged for the two 
separate fits II, and it might be tempting to interpret this as the O(a^) correction being 
very small and almost negligible. However, the central value for / in the combined fit 
is almost 10% smaller. 

We do not claim that the correct continuum limit can only be obtained with a combined 
fit. However, a proper continuum extrapolation of the results from the separate fits 
is not possible yet since data at two lattice spacings only are available. Moreover, 
the two lattice spacings do not cover a wide range with the fine lattice spacing being 
only 9 percent smaller than the coarse one. More data at significantly smaller lattice 
spacings seem necessary to shed more light on this issue. 

In section IIIIDI we argued that the F V corrections may be significantly larger for neutral 
pion masses much smaller than the charged ones. Having performed fits to the data we can 
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0.5798 


0.4682 


2.4901 


3.2702 


-0.462 


-0.150 


0.0035 


0.6999 


0.7722 


3.2370 


3.8693 


-0.157 


-0.070 


0.0055 


0.8090 


0.9369 


4.3628 


4.8504 


-0.040 
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0.8599 


1.0018 


5.2525 


5.6641 


-0.016 


-0.010 


0.0085 


0.8764 


1.0103 


4.2338 


4.5224 


-0.076 


-0.058 



TABLE IV: Results for the ratios Hogsi 
the combined fit V. 



"{MqL) and rM±L. Based on the parameters obtained in 



a posteriori quantify this enhancement. We define the ratio 

M2(ln(MoVAi)-^i(MoL) 

^logs 



Ml ln(M2/Ai)-^i(M±L) 



(59) 



of the neutral pion chiral log and the charged pion chiral log (with FV corrections included), 
and the ratio 



r(ML) 



~gi{ML) 
In(MVAi) 



(60) 



as a measure about the relative size of the FV correction to the (infinite volume) chiral log. 

Table |IV] gives the results for these ratios for the fit parameter -805^.3 and C2 of the 
combined fit V. Apparently, for the smallest quark mass (a/io = 0.0025) the neutral pion 
chiral log is less than half as large as the charged pion chiral log. And r(MoL) is three times 
larger than r{M±L). Although less severe we find a significant increase in the FV corrections 
for the heavier quark masses as well. 

The ETM collaboration has found in ref. j3] that the resummed FV corrections of 23 1 



describe the data better than the standard 1-loop FV corrections with the charged pion 
mass. This might not be an indication for the superiority of the resummed formulae. It may 
just signal the failure of the standard expressions due to the use of the heavier charged pion 
mass. 



V. CONCLUDING REMARKS 

As already mentioned, the numerical analysis of the previous section can certainly be 
refined in various ways. For example, the error analysis should be improved, and the proper 
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ratio of Z-f actors should be included instead of the approximate value 1 employed here. 
Nevertheless, the results of the previous section show strong evidence for the sensitivity of 
the data on the pion mass splitting in the chiral logs. The fits are in general better with the 
mass splitting taken into account. Moreover, provided the general trends of table UTTl persist, 
there is some sizable and non-vanishing impact of our new formulae on the extraction of 
LECs, in particular for / and I3. We also found that separate fits at fixed lattice spacings 
can be quite misleading concerning the scaling violations in these LECs. These findings, if 
corroborated, are of course very important for ChPT phenomenology. 

Further improvement seems possible if data for the neutral pion mass are included in 
the simultaneous fit. Although the neutral pion mass data are in general afflicted with 
significantly larger errors, including the data may still reduce the error on C2 together with 
improved estimates for the other fit parameters, in particular the physical LECs. 

An immediate question is whether more observables are affected in an analogous way 
by a large pion mass splitting. Other purely pionic observables are pion scattering lengths. 
These have been studied in twisted mass WChPT 45|], but only in the GSM regime where 
the pion mass splitting is a NLO effect and therefore ignored in the pion loops. 

More relevant in practice is the extension of the results given here to WChPT including 
the kaons. After all, the simulations in ref. 0] take into account a dynamical strange quark 
and we expect similar modifications in the 1-loop results for the mass and decay constant 
of the kaon. 
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